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Abstract. The effect of Rashba spin-orbit interaction in quantum wires with hard- 
wall boundaries is discussed. The exact wave function and eigenvalue equation are 
worked out pointing out the mixing between the spin and spatial parts. The spectral 
properties are also studied within the perturbation theory with respect to the strength 
of the spin-orbit interaction and diagonalization procedure. A comparison is done 
with the results of a simple model, the two-band model, that takes account only of the 
first two sub-bands of the wire. Finally, the transport properties within the ballistic 
regime are analytically calculated for the two-band model and through a tight-binding 
Green function for the entire system. Single and double interfaces separating regions 
with different strengths of spin-orbit interaction arc analyzed injecting carriers into 
the first and the second sub-band. It is shown that in the case of a single interface the 
spin polarization in the Rashba region is different from zero, and in the case of two 
interfaces the spin polarization shows oscillations due to spin selective bound states. 
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1. Introduction 

The Spintronics [I] is one of the most prominent fields of modern condensed matter 
physics. Its target is to use spin to create electrical and optoelectronic devices with new 
functionalities [2]. Up to now, starting form the seminal device by Datta and Das [3], 
several devices based on the giant and tunnel magnetoresistance have been realized for 
read-head sensors and magnetic random-access memories [2]. The important task of 
the integration of such spintronics technologies with the classical semiconductor devices 
finds an obstacle in the small spin injection from magnetic to semiconductor materials 
due to the large resistivity mismatch between magnetic and semiconductor materials pE]. 
For this reason it is useful to design semiconductor devices with an efficient all-electrical 
spin-injection and detection via Ohmic contacts at the Fermi energy, as it has been 
already realized for metallic devices [6] . 

Two important classes of spin-orbit interaction (SOI) are relevant for semiconductor 
spintronics: the Dresselhaus type [7j and the Rashba type [8] coupling. The former 
arises from the lack of symmetry in the bulk inversion whereas the latter arises 
from the asymmetry along the growing-direction-axis of the confining quantum well 
electric potential that creates a two-dimensional electron gas (2DEG) on a narrow- 
gap semiconductor surface. Since the Rashba SOI can be tuned by an external gate 
electrode 0, [TO, ELL] it is envisaged as a tool to control the precession of the electron 
spin in the Datta-Das proposal for a field-effect spin transistor [3]. 

Quasi-one-dimensional electron gases or quantum wires (QWs) are realized by 
applying split gates on top of a 2DEG in a semiconductor heterostructure [12]. The 
main effect owing to the confining potential is quantization of the electron motion 
in the direction orthogonal to the wire axis. The combination of this confining 
potential and the Rashba SOI gives rise to sub-band hybridization that can affect the 
working principle of the field-effect spin transistor. Mireles and Kirczenow p3] have 
numerically studied this effect and they have shown that a large value of the Rashba 
SOI can produce dramatic changes in the transport properties of the device till to 
suppress the expected spin modulation. The effect of sub-band hybridization has been 
investigated by Governale and Ziielike [H] in a QW with parabolic confinement. They 
show that electrons with large wave vectors in the lowest spin-spit sub-bands have 
essentially parallel spin. But in proximity of the anti-crossing points due to the sub- 
band hybridization it is no more appropriate to use the spin quantum number in order 
to characterize the electron state in the QW. Furthermore, they show that it is not 
possible to transfer the finite spin polarization of the QW to some external leads. 

In this Article we study the spectral and the transport properties of a QW in 
the presence of SOI. The main result of this Article is to show how to achieve non- 
zero spin polarization in external leads using spin unpolarized injected carriers. This 
can be obtained by injecting carriers in all the active sub-bands due to the quantum 
confinement. At the opening of each new sub-band the hybridization owing to SOI gives 
rise to spin selective bound states reflecting in a oscillating spin polarization. Here, we 
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want to stress that those spin polarized bound states are not in contradiction with any 
fundamental symmetry property of the system [15] . 

This Article is organized in the following way: in Sec. [2] we evaluate the spectral 
properties using the wave function approach [16j H3 EH]- m Sec. |3] we provide an 
exact calculation for the spectral properties investigated within the perturbation theory 
approach and with the exact diagonalization in a truncated Hilbert space. Here we also 
introduce a minimal model featuring the basic characteristic of a QW with SOI named 
two-band model. This is used in Sec. H]in order to study the transport properties of a 
QW in presence of a single interface between a region with and without SOI and in the 
case of a double interface (spin-field effect transistor scheme). Conclusions are ending 
the Article. 



2. Exact solution of Rashba quantum wire: wave-function 

Let us consider a 2DEG filling the plane (x, z). The charge carriers have momentum 
p = (px,Pz) and effective mass m. The particles are confined along the z-direction 
by the potential V(z) and subjected to the Rashba spin-orbit interaction (SOI). The 
single-particle Hamiltonian reads 

H = ^(pl+p 2 z )+V(z)+H R , (1) 
where is the Rashba SOI 

ri R = {a z p x - a x p z ) . (2) 

m 

In Eq.(J2J) o~ x and a z are the x and z components, respectively, of the vector a of Pauli 
matrices, and k$o is the SOI constant. This can be tuned by means of external gates 
perpendicular to the 2DEG P QUI [II]. 

In the following we assume that the potential V(z) provides a confinement with hard 
walls at z = and z = W. The strategy to find the wave- function of the Hamiltonian 
(JT|) is similar to the procedure followed in the absence of SOI: one exactly solves the 
2D problem, then considers the quantizing effect of confinement on the wave function 
ijj(x, z). In the first subsection we shortly recall the spin-dependent solution of the 2DEG 
with SOI, then we impose the boundary conditions ip(x, z = 0) = ?p(x, z = W) = 0. In 
the presence of SOI, these relations mix the z part of the wave-function with its spinor 
component. 



2.1. Solution without confinement 

Without confinement both components of the momentum p = hk are conserved. The 
eigenfunctions of the Hamiltonian (OQ) in the absence of confining potential V(z) = are 
denoted by the two spin modes (+) and (— ) and read 

*Pk,+ ( x > z ) = ex P Kkcx + k z z)\ ( ^^2) ) ' ( 3fl ) 
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Figure 1. (Left Panel) Spectrum of the two-dimensional electron gas in the presence 
of SOI interaction as a function of the wave- vectors k x and k y . At fixed positive energy, 
the wave- vectors of the mode (+) and (— ) are on two concentric circles (radius for ( — ) 
mode larger than that for the (+) mode). (Right Panel) At fixed positive energy and x 
component of the wave-vector, the 4 possible values of the z component are shown. 



(x, z) = exp [i(k x x + k z z)\ ^ ^j^j j ' ^ 
whose corresponding eigenvalues are given by 

E± = ^ (k 2 ± 2kk so ) , (4) 



where k = \Zk% + k\ is the modulus of the wave- vector in x-z plane, k x = kcos(8), 
k z = ksin(6), with 9 the angle formed between the vector k and the x axis. It is clear 
for Eqs. (13 al!3 b\ that the spinors \± °f t ne two modes are orthogonal to each other. 
We remind that the Rashba SOI can be viewed as a magnetic field parallel to the (x, z) 
plane and orthogonal to the wave-vector k. The net effect is to orientate the spin along 
the direction perpendicular to the wave- vector [46]. 

In order to determine the solution with confinement, it is essential to find the 
eigenfunctions in the free case when the total energy E and the momentum along k x are 
fixed (See Fig. [[]). Fixing the total energy E , we note that there are two values of total 
momentum k, corresponding to the different modes, fulfilling the Eq. (TjQ) expressed as 
a linear combination of those four waves: 



2m 

k± = \J-^E + kl Th . (5) 

The propagation directions for the k± modes are fixed by the momentum k x . For 
E > 0, the mode (+) is characterized by the propagation direction ±a = &rccos(k x /k + ) 
fixing the value of k z = ±A; + sin(a), whereas the mode (— ) has propagation direction 
±/? = arccos(fc x /p_) and k z = ±/c_sin(/3). Therefore, the generic wave function 
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^ E ,k x (x,z) is given by 



1> Bjhm (x, z) = e lfe ^ A^ E (z) + B^Q (z) + Cj>tt (z) + (z) , (6) 



with 



„/.(+) ( z) _ _-i(-l)<fc +2 sin(a) ( COs(a/2) \ 

^> e[ j_e ^ (-l)'siii(a/2) J' 10 

^C-)^ _ e -H-i)'h-z**w ( -(-l)'sin(/3/2) \ 

and £ = 1,2. 

For —k^<k x <—k + oik + <k x <k-, the wave-function ([6]) is still valid. However, 
one has a = ia, implying that cos(a) = cosh(a) and sin (a) = isinh(a), therefore the 
(+) mode becomes an evanescent one. Moreover, if k x > k_ or k x < —k_, then (3 = ib 
and also the (— ) mode changes into an evanescent one. 

For E < the 4 values of k z are only relative to modes (— ). A wave- function similar 
to ([6]) can be written. Also in this case the Fermi surface is formed by two circles but 
now they correspond to the same energy E_. However, in the next section, we will see 
that, from weak to intermediate values of the Rashba SO coupling, only positive values 
of the energy are important due to the effect of the confinement. 

2.2. Solution with confinement 

The wave-function (JHJ) represents the starting point for taking account of the 
confinement. In fact, the hard wall boundary conditions are obtained by imposing 
that the wave-function is zero on the borders {z — and z = W): ipE,k x {x,z = 0) = 
^E,k x {. x i z = W) = 0. For E > 0, we get the following exact eigenvalue equation for the 
Rashba quantum wire 

1 - cos[k + W sin(a)] cos[k-W sin(/?)] + 

sm[k + W sin(a)] sin[fc_Wsin(/?)] [1 + = 0- (9) 

sm(a) sm(p) 

Therefore, via the SOI, the quantities k + and k_, and clearly the energy, are related to 
the spinor components of the wave-function. A similar equation is valid for E < 0. 

In the absence of SOI, the quantized energy levels are independent of spin behavior. 
Actually, one gets k± = *J 2mE / h 2 and a = (3. This yields 



*2mE 

— — W sin (a) = mr, (10) 

h 2 

with n being a positive integer number, so that, together with the relation cos(a) = 
k x /k±, we obtain 

which are the energy values for the sub-bands of the quantum wire without SOI. 
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In the presence of SOI, an important limit is obtained for k x = 0. Indeed the 
eigenvalue equation becomes 

cos[(A; + + k-)W] = 1, (12) 

implying that 

Therefore, all the sub-bands are shifted down by the SOI term k^ Q . For values of the 
SOI such that kso < n/W, the energies are positives and the wave function (jfJJ) holds 
true. This means that the spin-precession length Lso = tt/^so has to be larger than 
the wire width. 



3. Two-band model and perturbation theory 

In the previous section we started from the wave-function of the 2D model with SOI, and 
then imposed the conditions due to confinement. Now, we consider the opposite point of 
view. First we take into account the exact solution of the quantum wire in the absence 
of the SOI, then we study its effect on the sub-bands. Because of SOI, a coupling 
between sub-bands with opposite spins occurs. In order to study the effects of this 
coupling, in the first subsection we will discuss the results within the first- and second- 
order perturbation theory approach with respect to the SOI. In the second subsection 
we consider the two-band model, where only the first two bands of the unperturbed 
spectrum are assumed to be coupled by the interaction. This assumption is valid if the 
wire is very narrow. Moreover, this simple system is studied since it provides a simple 
understanding of the transport properties. Finally, in the third subsection, we will show 
the results of the exact diagonalization of the model. 

The Hamiltonian (PQ) is considered to be split into two terms: 7Yo and 7Yr. The 
term Tio is simply the Hamiltonian of the wire without SOI: 

H = 7^{p 2 x + Pl)+V(z), (14) 

where V(z) is the hard-wall confining potential. Due to the presence of the potential 
V(z), only the momentum p x = hk x is conserved. We find the matrix elements of 7i 
in the basis of Tio indicated by \k x ,n,a), with n index of the sub-band and a = ±1 
for up or down spin, respectively. The SOI term contains terms o~ z p x and —o~ x p z . For 
£ = n, only the former term of 7Yr is acting on the unperturbed states with the same 
spin state, so that the matrix elements are 

(k x J,a\n R \k x ,n,a') = - ks ° kx a'8 (15) 

m 

while, for I ^ n, the latter term 7i R couples sub-bands with opposite spin and parity, 
so that the matrix elements are 

(k x ,£,a\H R \k x ,n,a') = J e>n 6 a _ al , (16) 
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with Ji n independent of the wave-vector k x 

= [i - (-i)"-"'] • (") 

m W £ 2 - n 2 1 J 
If we express the energies in the unit h 2 /2mW 2 , the lengths in W and the wave- vectors 
in 1/W, we recast the following matrix elements for the entire Hamiltonian "H: 

(k x ,£,a\H\k x ,n,a') = [E n °\k x ) + 2k so k x a'] 5<> tn 5 a y + 

Jl n [l - 8£ in ]5 a - a >, (18) 

where k x = k x W, E n °\k x ) = k 2 + n 2 ir 2 , kso = ksoW, and J^ n proportional to the 
dimensionless SOI term fcso 

A fin 

J*,n = ^SO^ [1 - (-1)1^1] • (19) 

3.1. Perturbation theory 

The correction to the unperturbed energies within the first-order perturbation 
theory is simply derived considering only the diagonal terms of Eq.( fl8l) . Therefore, 
at first-order, the n-th sub-band is simply affected by the spin slitting due to the 
contribution a z p x of Hso'- 

E n l l{k x ) = E n °\k x ) + 2k so k x a, (20) 

and eigenvectors equal to those of the unperturbed system. This splitting controlled by 
SOI gives rise to a first-order spectrum with crossings between sub-bands with opposite 
spins. For example, the first and second sub-band intersect at k x = ±37r 2 /4/cgo and 
the others for larger values of k x . This suggests that the full effect of the interaction 
should remove this crossing by mixing the behavior of coupled sub-bands. Due to the 
presence of those level crossings, the correction to the energy levels within the second- 
order perturbation theory fails for values of k x close to intersections. Far from the 
crossing points, it is easy to derive the contribution in the second-order to the energy 

^nl = Yl Yl K 2 \n 2 ^- I 2 ) = ^ 

a quantity independent of k x , n and a. This result is indubitably valid for k x = 0. 
Indeed, it coincides with the result (fT3l obtained in the previous section by using the 
exact wave-function. This shows that at k x = the energy correction within the second- 
order perturbation theory is able to fully describe the energy spectrum. Also, the 
correction of the wave-function at first-order can be evaluated. If at zero-order the spin 
is a, at first-order one takes contribution from —a: 

= i5Vw [SiAz) - S2Az)] 1 _ a) ' (22) 

with 

S 1>n (z) = $ ( e ~ 2mz /\ 2, ^ )- $ fe" 2 ^/ L , 2, , (23) 
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and 

S 2 ,n{z) = $ [e?^l\ 2, ^\ - $ (e 2 "^, 2, |) , (24) 
where s, a) = J^fcLo ( a + * s ^ ne Lerch transcendent function [19J. 
3.2. Two-band model 

In order to investigate the effects of the coupling between sub-bands induced by SOI, it 
is convenient to analyze the two-band model that will be also considered in the section 
devoted to transport properties. This model takes only the first and the second sub- 
band of the unperturbed wire into account . The 4x4 problem can be decoupled 
into two 2x2 problems. The only thing to evaluate is Ji t2 = — il6&so/3. One gets 4 
eigenvalues [T3| 



with 



e 1+ (k a 
9i(k x ) 



5tt2 

~Y 

5tt 2 



+ k\ - gi(k x ), €i_{k x ) = e 1+ (-k x ), 



k 2 x - g 2 (k x ), e 2 -{k x ) = e 2+ (-k x ), 



Ak^k 



, 1024^2 .- 

>r H g J 92\K 



<7i( 



(25) 
(26) 

(27) 



The eigenvectors can also be calculated. For example, the eigenvector corresponding to 
e 1+ is 



ip 1+ (x,z) 
with fi (k x ) given by 

h (h) = 




1 



lQk 



so 



w 



37T 5 



[A(k x )\ 



^ (f) 

if, {k x ) sin (^) 



(2? 



2k S ok x + gi(k x ) 



(29) 



As shown in Fig. |2j the eigenvalues (solid lines) do not show any intersection for 
k x different from zero. Therefore, the inter-band coupling removes the crossings of the 
first-order perturbation theory solution (dashed line). As a result, the energy eigenstates 
are no longer eigenstates of a z [H] and the spin state depends on the wave- vector k x . 
Close to the crossing point, the wave function of 1 f and 2 |, for example, are strongly 
mixed in the mode 1+. However, far from the intersection, the mode given by the 
diagonalization preserves the original behavior of the component wave-functions. For 
example, if we analyze the behavior of the eigenstate ipi + (x, z), we get 

lim i(; 1+ (x,z) = ip^(x,z), lim if) 1+ (x t z) = foifa z). (30) 

k x — >— oo k x —>oo 

The behavior of the two-band model shows a general trend: only taking into account 
the coupling between sub-bands the description is qualitatively correct. The crossing 
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Figure 2. Energy levels of the two-band model (in units of ti 2 /2mW 2 ) as a function 
of the dimensionless wave- vector k x W for the value k$oW = it of the dimensionless 
spin-orbit parameter. The spectrum derived from the diagonalization of the two-band 
system, indicated in figure by 1+, 1 — , 2+, and 2— (solid lines), is compared with that 
from first-order perturbation theory, indicated by 1 up, 1 down, 2 up, 2 down (dashed 
lines). 



are artifacts of the lowest-order perturbation theory. The spectrum within the two- 
band model is reliable only for very narrow wires. In the general case, the low-energy 
description given by this model is too poor for the bands 2+ and 2—. This can be easily 
seen if one considers the energy values at k x = 0. In fact one gets for sub-bands 1± and 
2±, respectively, the achieved values are 



5tt 2 1 / , 1024H 
ei± = — - 2 V 97r4 + — g-^ > ( 31 «) 

5tt 2 1 / ' 1024H 

In the limit of small k^o, they become 

ei± = " 2^io - ^ ~ 0.961^ o , (32a) 
256 

e 2 ± = 4tt 2 + i^klo - 4vr 2 + 0.961fc| o . (326) 

From the comparison with the exact solution we find that the lowest sub-bands acquire 
a correction with the right sign and very close to the exact result, while the upper sub- 
bands have even the wrong sign. Therefore, in order to give a reasonable description of 
the low-energy part of the spectrum, more bands are necessary. This will also play an 
important role in the transport properties. 



3.3. Exact diagonalization 



In order to verify the importance of including more than two sub-bands, one can directly 
diagonalize the Hamiltonian of the system [14] . This can be done considering the 
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-2% -71 71 2% 
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X 

Figure 3. Energy levels of the wire (in units of H 2 /2mW 2 ) as a function of the 
dimcnsionlcss wave- vector k x W for the value k$oW = n of the dimcnsionlcss spin- 
orbit parameter. The spectra derived from the diagonalization of the system with 3 
sub-bands (solid line) and 2 sub-bands (dashed line) are shown. 



matrix elements ([18]) . In Fig. [3] we consider the diagonalization in the subspace of 
3 spin degenerate bands. It is apparent that already at this level the corrections to the 
level energies and wave-functions are important for the second sub-band. For example, 
it is very close to the correct behavior at k x — 0. Considering I sub-bands for the 
diagonalization, one is able to get a reliable behavior starting from sub-band 1 to I — 1. 

4. Ballistic Transport 

In this section the issue is to study the quantum transport properties within the ballistic 
regime. The standard Landauer-Biittiker formalism will be employed. We will start 
considering a wire divided in two different regions: one with SOI and one without. 
Then we take the case of a quantum wire with a finite SOI region into account. The 
results are obtained within the approximation of the two-band model and those results 
will be compared with a numerical tight-binding method. 

4-1- Single interface 

We consider a QW divided in two main regions: in the right region (x > 0) SOI is present 
while in the left region (x < 0) it is not. The interface separating the two regions is 
considered to be sharp and is described by a 5-like potential. The Hamiltonian of this 
hybrid system reads 

^•hyb = P r, , s P + V[Z) H [pzPx ~ 0- x p z ) 

2m(x) m so 

. h dk so (x) h 2 u 
- mz 2m so dx (33) 
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We assume that the mass and the strength of the SOI are piecewise constant with 
kso{ x ) — kso@{ x )- F° r simplicity the mass is considered equal on both sides of the 
interface. The fourth term is necessary to get 7i ny b hermitian. At the interface the spinor 
eigenstates of 7Y ny b have to be continuous, whereas their derivatives have a discontinuity 
u — icr z kso due to the SOI and to the 5-like potential in x = 0: 

^(0+) = ^(0") , (34a) 



dij){x) 



dip(x) 



x =o+ dx 



= (u-\a z k so )^). (346) 

x=Q- 



dx 

In order to study the effects of the sub-band hybridization on the transport 
properties, we start considering the injection of carriers only within the first sub-band. 
This can be achieved requiring that the second sub-band is behaving as an evanescent 
wave. In this context the Eqs. f!34all346j) are reduced to a set of two decoupled systems 
of four times four equations for the variables ri+, t± + , r 2 _, t 2 ~ and ri_, fi_, r 2+ , 
t 2+ respectively, where ti( 2 ),+(-) an d r i(2),+(-) are the transmission and the reflection 
amplitudes in the first (second) sub-band with spinor + (— ) respectively. The knowledge 
of the transmission and reflection amplitudes permits to evaluate the probability current 
and, as a consequence, the transmission probabilities for spin-up and spin-down carriers. 
Due to the presence of SOI only for x > 0, the probability current has two different 
forms given by 

1 J 3? {ip^pip} for x < , . 

3 ~ m I 3? [p + hk so (y x S)\ i)} for x > ^ 

where i/j is the wave function solution of the system of Eqs. (l34aH346)) . A direct evaluation 
of the transmission probabilities [20] results in the following expressions: 



If 



12 



-!— [k 1+ (k in ) + k so (a z ) 1+ ] , (36a) 

rC;, 



If 



in 

|2 



T { = ^ [h-ihn) + fcso(Oi-] > (366) 

Kin 

where k m is the injection momentum and ki + (ki n ) and ki-(k m ) are the momentum 
relative to k[ n for the two spin resolved sub-bands. The factors (c z )i+ and are 
the expectation values of a z on the two spin resolved sub-band wave functions and are 
defined as 

{ ^ = ± i +M ± kl± r (37) 

where the function fi has been provided with Eq. (1291) . Because of the absence of 
SOI for x < the reflection probabilities are simply defined as = |r 1+ | 2 and 
R± = |r!_| 2 . The system Hamiltonian 7i ny b is invariant under time-reversal symmetry 
and, as consequence, the following relations hold: R^ = Ri and T-f = Tj_. This means 
that spin-up and spin-down carriers are transmitted through the interface in the same 
way. 
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Figure 4. Polarization as a function of the injection energy for fcgo = (dotted 
line), fcso = 1 (solid line), fcso = 2 (dashed line) and fcso = 3 (dotted-dashed line). 



We consider the injection of a spin-unpolarized mixture of carriers with injection 
momentum k- m . In terms of density matrix we have 

Pin = ^|T)(T| + ^|I)(!| (38) 

with the property that (cr z )i n = Tr{pi n a z } = 0. The density matrix of the transmitted 
carriers is expressed by the relation 



T 



Pout 



i+)(i + | + 



l+)(l + l + -ll- 



T t +T 



[1— ><1 



(39) 



\2 



(40) 



T " 1 2 

where |1±) are the wave functions of the spin resolved sub-bands. We can now evaluate 
the polarization of the output carriers, this is expressed by 

l i v 1 

i n-f 1 (k 1+ ) 2 l-M-h 
" 2 yi + M^+y i + h{-k^)\ 

Figure (HI) shows the polarization ( 14*01 as a function of the injection energy for various 
values of the SOI. The injection energy is limited within the first two spin-resolved 
sub-bands. It is clear that when SOI is zero (dotted line) there is no polarization, but 
as soon as the SOI is different from zero, the polarization gets a finite value, which 
increases as a function of SOI for a fixed energy. Those results are in accordance with 
Governale and Zulicke p3], who showed a negative polarization for carriers in the first 
two spin-resolved sub-bands and positive injection energies. 



4-2. Wire with finite SOI region 

We consider a QW composed of three parts: two external regions (x < and x > L) 
without SOI, and a central one (0 < x < L) where SOI is present. We consider the 
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hybrid system Hamiltonian of Eq. (1331) . with kso(x) = ksoO(x)9(L — x) and two 5-like 
potentials of x = and x = L. Also in this case, for simplicity, the mass is assumed 
constant along the wire. As for the single interface case, the spinor eigenstates of i^hyb 
are continuous at the interfaces, whereas their derivatives have discontinuities in x = 
and x = L: 



^(0 + ) 
difj(x) 



dx 
dip(x) 



dx 



^(o- 

=0+ 
=L+ 



dip(x) 



dx 



dx 



(u + ia z k S o)ip(L). 



(41a) 
(416) 
(41c) 
(Aid) 



As first step, we consider carriers with injection energy within the first two spin-resolved 
sub-bands and with evanescent waves for the following two. The Eqs. ( l41aH41t| reduce 
to a set of two decoupled system of equations, with relevant terms £ri+, £ri-, tli+ and 
rLi_. Those are practical for evaluating the transmission and the reflection probabilities 
for spin-up and spin-down carriers. As expected, because of the absence of SOI in the 
external regions we get 



1+1 



Ri 



\tm- 
-- ku- 



(42 a) 
(426) 



R\ = I'L1+| , -n-i 

Due to time-reversal symmetry, it results that the value of the transmission probability 
T| for spin-up incoming carriers is equal to Tj_ for incoming spin-down carriers. It is 
relevant to study transport properties when an unpolarized mixture of spin-up and 
spin-down carriers is injected into the system. The incoming and the outgoing density 
matrices are described by the expressions ( |38|) and ( l39l) respectively, where now the 
states |1±) are the output wave functions in the second region without SOI. As in the 
previous section we can evaluate the polarization as the average value of the a z operator 
and obtain as result 



(^)out = T| — T| = 0. (43) 

The effect of spin polarization due to the first interface is completely cancelled by the 
second one, therefore it is not possible to observe any spin polarization p3]. This result 
can be also derived by a symmetry consideration: let us consider the scattering matrix 
Sqw of the QW. In absence of a magnetic field time-reversal symmetry is preserved, 
therefore, as a consequence, for <Sqw holds the following important relation: 



where S, 



2 



0c 



0~v 



It is clear for Eq. ( 1441) that in the case of only one 



conducting channel the spin-flip transmission terms must be zero and as consequence 
the polarization is absent [T5] . 



Rashba quantum wire: exact solution and ballistic transport 

3, 1 1 1 1 1 1 1 1 01 1 ' — 



11 




f— 1- 



-0.5 



2.5 



Injection energy 



5 7.5 
Injection energy 



10 



Figure 5. (Left Panel) Spin resolved transmissions as a function of the injection 
energy and for two different values of the transparency of the barriers: it = 1.0 (solid 
lines) and it = 0.1 (dashed lines). The dotted line indicates the opening of the sub- 
bands in absence of SOI. The dimensionless SOI is fcso = 1.4 and the central region 
is L = 3W . (Right Panel) Spin polarization (a z ) as a function of the injection energy 
and for two different value of the 5-barrier transparency. The same parameters as in 
the left panel are used. 



As a further step, we study the transport properties when carriers are injected also 
within the second two spin-resolved sub-bands. In this limit all the evanescent modes 
are transformed in conducting ones. Using Eqs. (}38l) and (|39|) for the incoming and the 
out-coming density matrix, the polarization is expressed by 

_ 1 (T 1+ (0) + T 2+ (0) - Ti_(0) - T 2 _(0) 



+ 



+ 



T 1+ (0) + T 2+ (0) + Ti_(0) + T 2 _(0) 
T 1+ (7r/2) + T 2+ (7r/2) - ri_(7r/2) - T 2 _(tt/2) 



(45) 



T 1+ (tt/2) + T 2+ (7r/2) + T 1 _(7r/2) + T 2 _(7r/2). 
where "0" and "7r/2" denote incoming up- and down-carriers respectively. 

In Fig. [5] (left panel) we show the spin resolved transmissions as a function of the 
injection energy below and above the bottom of the second two sub-bands and for two 
values of the transparency u of the interfaces (strength of the ^-function potential in 
Eq. l33j) . The behavior below the threshold clearly shows Fabry- Perot oscillations due 
to multiple reflection effects within the central region, moreover the strength of those 
oscillations is decreasing for increasing transparency of the two 5-barriers. A second 
important feature of the behavior is that up to the threshold, owing to the time-reversal 
symmetry, the two spin resolved components possess the same value. When the injection 
energy is crossing the bottom of the second two sub-bands, a new behavior appears. The 
two spin resolved transmissions take different values, the difference between those two 
is bigger near to the sub-band threshold and is decreasing with increasing energy. This 
can be well understood if we relate this phenomenon to the sub-band hybridization. 
As shown in Figs. [2] and [3] the SOI strongly modifies the parabolic-like behavior of the 
spin resolved sub-bands when new sub-bands are opening and this effect is vanishing at 
higher energies. In contrast to the results of Governale and Zulicke [Hj, Fig. [5] clearly 
shows how polarization effects can manifest themselves when also the second two spin- 
resolved sub-bands are involved into the transmission mechanism. This is more clear 
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Figure 6. (Left Panel) Numerical spin resolved transmissions as a function of 
the injection energy and for two different values of the switching region of SOI: 
L sr W~ x ~ 0.07 (solid lines) and L^W^ 1 ~ 0.7 (dashed lines). The dotted line 
indicates the opening of the sub-bands in absence of SOI. The dimcnsionless SOI 
is fcso = 1.4 and the central region is L = 3W + 2L sr . In the insets a) and b) are 
shown the magnification for the steps relative to the opening of the second and the 
third sub-band, respectively. (Right Panel) Spin polarization (a z ) as a function of the 
injection of energy and for the two different values of L^W^ 1 . The same parameters 
as in the left panel are used. 

in Fig. [5] (right panel), where the polarization as a function of the injection energy is 
shown. The polarization is zero below the bottom of the second two sub-bands and 
has oscillating behavior above it. This particular shape of the polarization can be 
understood also in terms of formation of spin-dependent bound states at energies closer 
to the opening of the second two sub-bands. Those bound-state oscillations are clearly 
visible in the polarization pattern and it is evident how they are enhanced when the 
(^-barrier transparency is decreased. 

The effects on spin transport due to higher sub-bands can be numerically verified. 
For this purpose we employ a tight-binding model of the Hamiltonian (CD), the spin- 
dependent scattering coefficients are obtained by projecting the corresponding spin- 
dependent Green function of the open system onto an appropriate set of asymptotic 
spinors defining incoming and outgoing channels. A real-space discretization of the 
Schrodinger equation in combination with a recursive algorithm for the computation of 
the corresponding Green function has been implemented [20]. This formalism allows 
a convenient treatment of different geometries as well as different sources of scattering 
within the same framework [21~] . 

The <5-barriers are introduced in the tight-binding calculation dividing the central 
region in three parts in order to modulate the switching region L sr for the SOI. A long 
switching region will correspond to a high transparency, and vice versa. In Fig. [6] (left 
panel) we show the spin-resolved transmission as a function of the injection energy and 
for two different lengths of the switching region. The system parameters are chosen in 
order to have three active sub-bands in the case of injection within the highest energy 
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Figure 7. Polarization as a function of the injection energy for the two-band (dashed 
line) and the iV-band (solid line) models, respectively. Panel a) first sub-band, Panel 
b) second sub-band. For both the Panels N = 50. 



allowed by the tight-binding approximation of the Hamiltonian ([I]). Injecting carriers 
within the first two spin-resolved sub-bands reproduces the result coming from time- 
reversal symmetry: spin-up and spin-down transmissions coincide. As soon as the second 
two sub-bands are opened we can observe a difference in their values. This difference 
tends to decrease for increasing energy but is, then, enhanced by the opening of the 
third two sub-bands [22]. This is clear from Fig. [6] (right panel), where we show the 
polarization as a function of the injection energy. Here it is also possible to observe how 
the strength of the polarization oscillations in strongly reduced changing the length of 
the switching region. 

The polarizations at the opening of the second two spin-resolved sub-bands in the 
case of the two-band model (Fig. [5]) and for the numerical tight-binding method (Fig. [6]) 
show an opposite value. This can be well understood if we consider the polarization 
within the first and the second sub-band for an infinite long wire in the two-band and 
in the iV-band models, respectively. In Fig. [7] we show the polarization for the first 
sub-band (upper Panel), it is evident that there is a small deviation of the two-band 
model (dashed lines) from the iV-band model. But this deviation is stronger for the 
second sub-band (lower Panel) resulting in completely opposite value at large energies. 
This well reveals how the two-band model can only qualitatively reproduce the exact 
calculation, but is failing in the quantitative estimation. 

5. Conclusions 

We have studied the properties of a QW in the presence of Rashba SOI. The spectral 
properties have been investigated both from the point of view of the exact wave 
functions and within the second-order perturbation theory approach. Furthermore, a 
numerical diagonalization procedure has been implemented in order to study the spectral 
properties for the system in presence of N sub-bands. We have used this last method 
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with two sub-bands within the so called two-band model. We have shown that this is a 
good description for the full properties of the first sub-band but not of the second one. 
This model has also been used in order to study the transport properties of a system with 
an interface between a region with and without SOI. We have shown that the interface 
is spin selective and that the crossing of unpolarized carriers through the interface can 
give rise to a non-zero polarization. We have studied spin transport also injecting 
carriers within the first and the second sub-band. We have shown that at the opening 
of each new channel the sub-band hybridization can give rise to spin selective bound 
states reflecting in a oscillating spin polarization. Those results have been confirmed by 
numerical calculation obtained within the tight-binding approximation. 



Acknowledgments 

We acknowledge I. Adagideli, M. Scheid and M. Strehl for useful discussions. D.B. 
wants to thank the hospitality at the Max- Planck- Inst it ut fur Physik komplexer 
Systeme (Dresden) where this paper was partially written and the financial support 
of Deutsche Forschungsgemeinschaft within the cooperative research center SFB 689 
"Spin phenomena in low dimensions" . 



References 

I. Zutfc, J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 76, 323 (2004). 
S.A. Wolf et at, Science 294, 1488 (2001). 
S. Datta and B. Das, Appl. Phys. Lett. 56, 665 (1990). 
G. Schmidt et ai, Phys. Rev. B 62, 2790 (2000). 

F. J. Jcdema, A.T. Fillip, and B.J. van Wess, Nature 410, 345 (2001). 
S.O. Valenzuela and M. Tinkham, Nature 442, 176 (2006). 

G. Dresselhaus, Phys. Rev. 100, 580 (1955). 

E. Rashba, Fiz. Tverd. Tela (Leningrad) 2, 1224 (1960), [Sov. Phys. Solid State 2, 1109 (I960)]. 
J. Nitta, T. Akazaki, H. Takayanagi, and T. Enoki, Phys. Rev. Lett. 78, 1335 (1997). 
Th. Schapers, J. Engels, T. Klocke, M. Hollfelder, and H. Liith, J. Appl. Phys. 83, 4324 (1998). 
D. Grundler, Phys. Rev. Lett. 84, 6074 (2000). 

Th. Schapers, J. Knobbe, and V. A. Guzenko, Phys. Rev. B 69, 235323 (2004) 

F. Mireles and G.Kirczcnow, Phys.Rev.B 64, 024426 (2001). 
M. Govcrnale, U. Zulicke, Phys. Rev. B 66, 073311 (2002). 
Feng Zhai, and H. Q. Xu, Phys. Rev. Lett. 94, 246601 (2005). 

V. M. Ramaglia, D. Bercioux, V. Cataudella, G.De Filippis, C.A. Perroni and F. Ventriglia, Eur. 

Phys. J. B 36, 365 (2003). 
V.M. Ramaglia, D. Bercioux, V. Cataudella, G.De Filippis, and C.A. Perroni, J. Phys.: Condens 

Matter 16, 9143 (2004). 
D. Bercioux, and V. M. Ramaglia, Superlattices Microstruct. 37, 337 (2005). 
A. Laurincikas and Ramunas Garunkstis The Lerch Zeta-Function (Kluwer Academic Pub). 
D. K. Ferry and S. M. Goodnick, Transport in Nanostructures (Cambridge University Press, 1997). 
D. Frustaglia, M. Hentschel, and K. Richter, Phys. Rev. B 69, 155327 (2004). 
Bound states owing to Fano resonances in quantum wire with Rashba spin-orbit interaction have 
been recently investigated also by Sanchez and Serra, Phys. Rev. B 74, 153313 (2006). 



